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Abstract 

We investigate the relaxation of braided magnetic loops in order to find out how the type of 
braiding via footpoint motions affects resultant heating of the loop. Two magnetic loops, braided 
in different ways, are used as initial conditions in resistive MHD simulations and their subsequent 
evolution is studied. The fields both undergo a resistive relaxation in which current sheets form 
and fragment and the system evolves towards a state of lower energy. In one case this relaxation 
is very efficient with current sheets filling the volume and homogeneous heating of the loop oc- 
curring. In the other case fewer current sheets develop, less magnetic energy is released in the 
process and a patchy heating of the loop results. The two cases, although very similar in their 
setup, can be distinguished by the mixing properties of the photospheric driver. The mixing can 
be measured by the topological entropy of the plasma flow, an observable quantity. 
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1 Introduction 

Coronal loops are enormously diverse in their nature, acting as building blocks of the corona, from 
bright points to active regions and flaring loops. As such, loops cover a huge range of lengths 
(1 - 1000 Mm) and it seems likely that several coronal heating mechanisms are responsible for 
heating loops to the observed range of temperatures (0.1- >10 MK). Explaining the observations 
remains a challenge and a number of questions are currently under debate. For example, can 
loops be broadly classified as isothermal or multi-thermal (e.g. Schmelz et al. 2009; Aschwanden 
& Boemer 2011)? Is heating impulsive or steady (e.g. Patsourakos & Klimchuk 2008; Tripathu 
et al. 2010; Warren et al. 2010)? A recent review of observations and modelling of coronal loops 
can be found in Reale (2010). Furthermore, a more general coronal heating review is given by 
Khmchuk (2006). 

One very promising loop-heating method, following the early ideas of Gold (1964) and Parker 

(1979, 1994), is magnetic braiding. Here photospheric motions acting on the loop footpoints act to 
twist and tangle the overlying field, increasing its magnetic energy. Eventually current layers (sin- 
gular or non-singular) may form in the field (e.g. Longcope & Sudan 1994; Galsgaard & Nordlund 
1996; Longbottom et al. 1998; Ng & Bhattacharjee, 1988; Craig & Sneyd 2005; Wihnot-Smitii 
et al. 2009). Magnetic reconnection will then enable a restructuring of the field as it relaxes to 
a lower energy state with plasma heating a natural consequence of the energy release. Magnetic 
braiding is also a possible explanation for the observation that coronal loops have approximately 
constant width (Khmchuk 2000; Lopez Fuentez et al. 2008): a braiding of field lines within a 
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Figure 1 : Cartoons showing the way in which the fields (left) and S ^ (right) can be built up by 
rotational footpoint motions acting on initially straight strands. 



loop prevents the expansion seen in simple potential or linear-force-free models (Lopez Fuentez 
et al. 2006). 

Although it is well-established that there is, in principle, sufficient energy in these photo- 
spheric motions for this mechanism to be plausible (e.g. Klimchuck 2006), whether the process 
is responsible for the coronal heating depends on many, often unkown, properties of the driver 
and the relaxation mechanism in the corona. One important question is how efficient the surface 
motions are at building up free energy in the magnetic field. Another unresolved issue is whether 
the relaxation mechanism in the corona can release this energy again. A theory often invoked to 
describe the energy release process is that of Taylor relaxation. This theory was initially devel- 
oped for a laboratory plasma device (Taylor 1974, 1986) but has also been applied to the solar 
case (Heyvaerts & Priest 1984, Dixon et al. 1989, Nandy et al. 2003, Kusano 2005, Hood et 
al. 2009). Under this hypothesis, the field relaxes to a particular linear force-free field (with the 
same global helicity, toroidal flux and boundary conditions as the initial field) so that not all the 
magnetic energy in excess of that of the potential field can be released. 

The aim of this paper is to show that the amount of energy that can be released in any non- 
ideal relaxation depends greatly on the topological properties of photospheric flow, which in turn 
determine the way in which the magnetic field lines making up the flux tube are braided (mixed 
and tangled together). To show this we consider the resistive evolution of two contrasting magnetic 
fields, both generated through rotational motions on the boundary. The first field has a sequence 
of rotational footpoint motions of alternating sense. The comparison case is also generated by 
rotational footpoint motions but these are all in the same sense. 

Our approach differs from studies (e.g. Gudiksen & Nordlund 2002; Peter et al. 2004; Bingert 
& Peter 2011) in which quasi-stationary processes of continuous driving and relaxation are mod- 
elled in that we first assume an ideal process where the photospheric driver braids a coronal loop 
to a certain level before a resistive evolution is allowed. This has the advantage that the more 
complex structures we aim to study can be built up without being quickly dissipated by high 
numerical resistivities. By contrast the present study involves a less sophisticated treatment of 
certain physical effects such as heat conduction and radiative losses, with results not yet suitable 
for forward modelling. As such we view the two approaches as being complementary. 

An outline of the work is as follows: In Section 2] we introduce the two magnetic fields {E^ and 
S^) whose MHD evolution will be studied throughout. Tools for measuring the level of braiding 
are discussed. In Section [3] we describe the numerical methods used for the simulations. Results 
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are presented in Section[4]and are broken down into two parts, firstly an examination of the basic 
properties of the resistive relaxation and secondly energetic considerations including estimates for 
coronal heating. We conclude in Section [5] 

2 Model Magnetic Loops 

Throughout this work, we employ an idealised representation of the coronal geometry, with nearly 
straight coronal fields running between two parallel places, which represent the photosphere. The 
first of our two magnetic fields, , is based on the pigtail braid. It may be built up from a uniform 
vertical field by rotational stirring motions on the boundary (the photosphere). The manner in 
which these motions act on uniform strands to create the braid is illustrated in Figure [l] (left). 
Two regions of rotational footpoint motion are present. As viewed from the direction shown in 
the figure, the right-hand motion acts in a clockwise direction and the left-hand motion in an 
anti-clockwise direction. Each motion rotates the strands about each other by a relative angle of 
n radians. The rotations are applied in the sequence o-y,^^^ ,o-\,o-2^ ,o-\,o-2^ , which is also the 
braid word (Birman 1975) representing this braid (the power negative one indicates the change 
of orientation of the rotation). Since the number of left and right hand rotations is the same the 
total magnetic helicity of the configuration is zero. While the cartoon image shows just three field 
lines of the braid, the motions on the boundary will affect all field lines within a certain range and 
generate a continuum of braid patterns. Those field lines lying outside the domain of the rotational 
motions will remain straight and undisturbed. 

Our second magnetic field, labelled 5 ^, is built up by a very similar sequence of motions on the 
boundary but with rotations all acting in the same, clockwise direction (cTj" \ cr" ^ , o"~ \ ct" ^ crj ' , ct" ' ). 
The braid representation of this field is shown in Figure [T] (right) where we see that the corre- 
sponding three strands have each undergone exactly 2n rotation. Again a continuum of braiding 
patterns will be found in all field lines making up the loop but the total helicity of the field no 
longer vanishes. 

We now construct an explicit magnetic field representation of these idealised pictures. As- 
suming a Cartesian geometry, we take a uniform background field (le^) and superimpose six flux 
rings, evenly spaced in z and located alternately at {x,y) = (1,0) or ix,y) = (-1,0). Each ring has 
components only in the Cv and Cy directions and is localised in all three dimensions. Together this 
creates six localised regions of twist in an otherwise uniform field. The closed form expression 
which generates such magnetic fields is given in Wilmot-Smith et al. (2009). For E^ we take 
exactly the same parameter set as given in that paper while for 5 ^ we simply change the sign of 
the twist parameter {k) to be the same {k - 1) for each twist region (rather than alternating, as for 

The aspect ratio of both loops is high (1 : 8). This is broadly consistent with observations 
of coronal loops whose elementary strands (width < 2 Mm, Aschwanden 2005) are much longer 
than they are wide. Both initial magnetic fields contain the same amount of magnetic energy 
( B^/2/j.o dV). The braiding has been applied in a conservative manner with the magnetic energy 
being only a small amount (3.08%) above that of the uniform background field. In other respects 
the two braids are very different, as detailed in the following paragraphs. 

To understand these differences consider the nature of the field line connectivities. An estab- 
lished method for doing so is to examine the squashing factor Q (Titov et al. 2002). This is shown 
on the lower boundary of both fields in Figure |2] While the maximum value of Q in both fields 
is comparable (specifically, Q{E^)max = 2.3 x 10^ and QiS^),nax - 2.4 x 10^), it is clear that 
there are many more layer-like regions of high Q for E^ than for 5^. These regions arise from 
the property that field lines making up E^ have a more complex connectivity than those of 5^, 
there being many more regions in which neighbouring field lines diverge as they are traced up 
through the corresponding braids. Overall this gives a 'mixing' of field lines with respect to their 
connectivities on the lower boundary. Simplistically, the efficiency of the field line mixing could 
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Figure 2: Distribution of logio(Q) lower boundary for the initial state of (left) and S ^ (right). 
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Figure 3: Colour maps in the initial states of (left) and (right). 

be quantified by, for example, integrating Q over the surface. Calculating (3 = jT logjQ(2)<iA 
(where A is the surface [-3, 3] x [-3, 3] shown in Figure |2]) for both fields we find Q^3 = 89.5 
while Qg3 = 68.0. These values confirm the qualitative picture given in Figure|2]that the field line 
mixing is better for E^. 

A second way to examine field line connectivity and mixing is by way of colour maps (Polymilis 
et al. 2003) as shown in Figure [3] . These will be particularly instructive later when considering 
how the systems relax to lower energy states. To explain how these images are generated, con- 
sider a field line threading the domain, and let the points of intersection of the field line with the 
lower and upper boundaries be (xo,3'o) and {X, Y), respectively. We make a plot over the lower 
boundary, colouring the point (xo,3'o) red ifX > xq and Y > yo, green if X < xq and Y < yo, blue 
if Z > xo and Y < yo and yellow if X < xq and Y > y^. In this way each point on the boundary is 
coloured, excepting those periodic orbits (X = xq, = yo) which (generically, as in these fields) 
lie at the intersection of all four colours. The complex colour maps for the initial states of 
and S ^ show the braided nature of the fields, E^ being more complex with small-scale structures 
filHng a greater portion of the domain. 

The colour map is a visual representation of the complexity in field line connectivity. A more 
precise, formal measure, is given by a quantity known as the topological entropy. Somewhat like 
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our integrated squashing factor, this is a global measure that gives a single number for the whole 
magnetic field. It has the advantages both of a firm theoretical grounding and of being a robust 
quantity insensitive to small changes in the magnetic field. There are several equivalent definitions 
of the topological entropy, but a convenient one is the asymptotic growth rate (with z) of horizontal 
loops stretched around the magnetic field lines (Newhouse & Pignataro 1993; Thiffeault 2010). 
While the exact entropy depends on the full pattern of magnetic field lines, a good estimate may 
be obtained by ensemble averaging over finite sets of field lines. Applying the numerical method 
of Moussafir (2006), as implemented by Thiffeault (2010), and with sets of 40 field lines, we find 
r(£'3) « 3.3 and T{S^)x23. 

In summary, while the two magnetic fields and S ^ may be generated by the same amount 
of boundary motion, the details of the pattern of boundary motions is crucial in determining the 
braiding pattern of the resulting fields. From the measures presented above we see that has a 
significantly higher degree of complexity than 5 ^*. With these differences in mind we now wish 
to use these two magnetic fields as initial conditions for resistive MHD relaxations. We aim to 
show that the ability to effectively heat large regions of the loop depends on the nature of the 
photospheric driver. Before presenting results of these experiments we first proceed to detail our 
numerical methods. 

3 Methods 

The magnetic fields corresponding to the closed-form expressions for E^ and S ^ are not force- 
free. However, the corona itself is thought to be a largely force-free environment, i.e. one with 
magnetic fields B and associated currents (J = (V x B)//io) satisfying J x B w 0. The first stage of 
our experiments is therefore to use an ideal Lagrangian relaxation code (Craig & Sneyd 1986) to 
relax the fields towards a force-free equilibrium whilst exactly preserving their topology. In the 
scheme, details of which may be found in Craig & Sneyd (1986), an artificial frictional evolution 
is taken to minimize J x B. The relevant output is then the final state of relaxation while the path 
to this state is not important. 

Each of the fields E^ and detailed in Section [His used as an initial condition for this ideal 
relaxation, over a domain x,y € [-6,6], z e [-24,24] with a uniformly spaced grid of 101^ 
points. The result of this procedure was described in detail for E^ in Wilmot-Smith et al. (2009). 
A smooth near-force-free field is obtained with large-scale current distributions in the form of two 
tubes of current running through the domain (see Figure [4]). By 'near-force-free' we mean that 
the maximum Lorentz force in the domain is (J x B),„av ~ 0.059 where both (9( J) and (9(B) ~ 1 . 
Numerical difficulties, as documented by Pontin et al. (2009) prevent possible relaxation to a field 
arbitrarily close to J x B = 0. For our purposes it is enough to be close to (rather than exactly 
at) a force-free state since this is also the relevant case in the solar corona. The ideal relaxation 
procedure for 5^ also results in a smooth approximately force-free field (with (JxB)maA: ~ 0.057). 
A large-scale current distribution is present, now with one twisted current tube running through the 
domain (see Figure |4]). In both cases the magnetic energy of the ideally relaxed field is reduced, 
the amount of energy in excess of potential now being 1.286% for E^ and 1.178% for . This 
will be discussed further in Section [421 

These approximately force-free fields are now used for the main body of our work in which 
they are taken as initial conditions in 3D resistive MHD simulations. In order to create the ini- 
tial conditions on the regular grid required, an interpolation procedure must be followed. The 
procedure is detailed in Wilmot-Smith et al. (2010) and ensures the interpolated field remains 
divergence-free to accuracies of the order of the truncation errors for sixth-order finite differences 
( |V • Blmax ~ 10"^ within the domain). We use the colour-map technique to check the conserva- 
tion of field line connectivity in this interpolation step and find the conservation to be good. 

The computational setup for our MHD simulations is described below. The experiments are 
conducted using the 3D non-ideal MHD code of Nordlund & Galsgaard (1997). This is a high 
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order finite difference code solving the following set of equations: 



-VxE, (1) 



dB 

~dt 

E - - (V X B) + 7/J, (2) 
J - VxB, (3) 



dp 
dt 

d 
dt 

'dt 



-V-(pv), (4) 
(pv) - -V • |pvv + rj - VP + J X B, (5) 
-V • (ev) - P V • V + Q,i,c + Qj, (6) 



where B is the magnetic field, E the electric field, v the plasma velocity, 77 the resistivity, J the 
electric current density, p the density, t the viscous stress tensor, P the pressure, e the internal 
energy, Qvisc the viscous dissipation and Qj the Joule dissipation. An ideal gas is assumed, and 
hence P = (y - \) e = |e. These equations have been made dimensionless by setting the 
magnetic permeability /uq = I, and the gas constant CR) equal to the mean molecular weight (M). 
Accordingly time units are such that, for a volume with |p| = |B| = 1, an Alfven wave would travel 
one space unit in one unit of time. 

We solve the equations over a grid with 256^ nodes over x,y € [-6, 6], z € [-24, 24], though 
during the simulations we find the dynamics to be confined approximately within x,y e [-4, 4] . 
The magnetic field is line-tied on all boundaries throughout, and the plasma velocity is fixed to 
zero at these boundaries. We obtain our initial magnetic field for the simulations via the interpo- 
lation method described above. The dimensionless plasma density is initialised as p = 1 and the 
thermal energy as e = 0.1. A spatially uniform resistivity model is taken, with t] = 0.001 for both 
simulations. 



In the energetic considerations of Section 4.2 dimensional quantities are recovered in order 



to demonstrate the implications of our results for the solar corona. In order to do this three 
characteristic values (here Bq, Iq and po) should be chosen and the following relations taken: 

vo = -, Bo^ vo V)"oPo, Jq = £'0 ^ VqBq, 
to fJ^oh 

eo = PqVq, Tq = —, 
K 

with = 47rx 10~^ H m"^ /i = 0.6 and /? = 8.3 x 10^ m^ s"^ (we use mks units throughout). 
For clarity we initially present results in the dimensionless units to allow the reader to adjust the 
chosen solar parameters as desired. 

Isosurfaces of current in the initial states for the resistive MHD simulations of and are 
shown in Figure[4] The current isosurfaces are at |J| = 0.441 for and |J| = 0.489 for 5^. These 
values may be compared with typical magnetic field strengths of |B| - 1 showing that, in both 
cases, these current ribbons are weak. 



4 Results 

The basic properties of the resistive relaxation of E^ have been described in Wilmot-Smith et 
al. (2010) and Pontin et al. (201 1). The braided field is found to be unstable (although the precise 
nature of the instability is yet to be determined) and in the early stages of the resistive evolution 
two main thin current layers are formed in the central regions of the domain. These fragment into 
an increasingly complex pattern of current layers as the relaxation proceeds. The lower energy 
end-state of the relaxation is found to consist of two unlinked twisted flux tubes of opposite twist. 
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Figure 4: Isosurfaces of | J| in the initial state of the resistive MHD simulations for (left) and S ^ 
(right). The isosurfaces are taken at 25% of the domain maximum in both cases. 



A resistive relaxation is also found to take place for the field . In this paper we aim to 
compare the two relaxation events in order to determine how the initial field configuration affects 
the quantity and spatial distribution of energy released and the final equilibrium state reached. We 
begin by discussing the basic nature of the relaxation events and final states (Section \A.\) before 



discussing findings related to energy and heating in Section 4.2 



4.1 Topological properties 

In qualitative terms the resistive MHD evolutions of and 5 ^ are somewhat similar. The large- 
scale currents of the initial state (Figure [i]) collapse towards two thin current layers for and one 
thin current layer for 5^. This process is shown in Figures [5j|6] and [v] It is at this early stage that 
the strongest currents are present in the system (see Figure [8j upper panel). As time progresses 
more current layers are found which are, individually, weaker in their intensity. For E^ a large 
number of current layers appear and they have a volume-filling effect (see, for example, f = 50 in 
the relevant figures). For 5^, by contrast, a smaller number of current layers are present and these 
are mainly located in the central regions of the domain (with respect to x and y). To quantify this 
diff'erence we show in Figure [s] (lower panel) the fraction of the volume (in [-4, 4]^ x [-24, 24]) 
for which |J| > 0.5 for both E^ (solid line) and (dashed line). The quantity is consistently 
larger for E^ , reaching a maximum of 8.1%, compared with only 2.8% for . In the later stages 
of the evolution the systems move towards equilibria containing only large-scale, weak currents. 
In the case of E^ these form two vertical tubes of current while for S ^ only one current tube is 
present. Note that we have run these simulations for a long time (up tot - 650) and find the state 
reached at f = 350 represents the end-state since after this time only slow changes in the magnetic 
field take place and these are due to global diffusion. 

Considering the final states of the resistive relaxations, although the number of current tubes 
present in each case is the same as in the respective initial states, the structure of the magnetic 
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Figure 5: Isosurfaces of current at 50% of the domain maximum for (upper panel) and (lower 
panel) at times t = 10, 50, 80, 120 and 350. 
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Figure 6: Contours of the vertical (e,) component of current in the z = plane for at times as 
indicated in each individual image. 
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Figure 8: Upper panel: Maximum value of the current density |J| in the domain with time. Lower 
panel: Fraction of the volume over which |J| > 0.5. In both cases the solid line corresponds to and 
the dashed line to 5^. 




Figure 9: Some illustrative field lines in the final states for (left) and (right). The field lines 
have been chosen to show the fundamental structure of the states, two flux tubes of opposite twist for 
E^ and a single twisted flux tube for 5^. 
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Figure 10: Colour maps in the final states of the resistive MHD simulations for (left) and 5^ (right). 



fields is very different. The end state for consists of two unlinked magnetic flux tubes of 
opposite sign of twist, as shown in the left-hand image of Figure [9] The end state for , by 
contrast, consists of a single magnetic flux tube of positive twist as shown in the right-hand image 
of Figure |9] The particular field lines plotted have been chosen because they portray the nature of 
the corresponding continuous magnetic fields. We show this nature in two ways, firstly in Figure 
10 using the colour-map technique described in Section|2] The corresponding colour maps for the 



initial states of both fields are shown in Figure [3] The first, striking, feature seen is the simplicity 
of the colour-maps in the final states in comparison to the initial states. This simplification comes 
from the un-braiding of the coronal loops to form two flux tubes for and one for S^. The 
difference in number of flux tubes may also be seen in the colour maps, there being two distinct 
intersections of all four colours (periodic orbits) for E^ and only one for 5^. These intersections 
mark the centre of the flux tube axes. 

A second feature of the colour map tells us why we cannot relax to a single flux tube in the 
E^ simulation (as determined by Yeates et al. 2010 , Yeates & Hornig 201 1). If we walk ai^ound 
the boundary of the domain for E^ then we meet each of the four colours twice (and in the order 
red-yellow-green-blue). This is because the overall topological degree of the state is +2. This 
is the same boundary pattern that was present in the initial state of E^ since the relaxation only 
affects the interior of the domain; mathematically, the topological degree of the field is conserved 
in the relaxation. Furthermore, a topological degree of +2 is inconsistent with a single flux tube. 
Considering S^, the boundary of the domain shows each colour appearing only once (but in the 
same sequence, topological degree -i-l) and so relaxation to a single flux tube is possible. 

A second way to examine the nature of the final state is to look at the mean value along field 
lines, a* , of the quantity a* - J • B/B • B. Figure 11 shows a* on the lower boundary of the 
domain for both the initial and final states of the resistive MHD simulations. For a perfectly 
force-free field, a* is simply the force-free parameter a. In our simulations the initial state is only 
approximately force-free and the finite gas pressure gives relaxed fields that are also not perfectly 
force-free. In both cases a* is the appropriate quantity to consider. The images in the initial state 
further illustrate the complexity of the magnetic field while the images for the final state confirm 
them as being much simpler. The basic un-braiding into two {E^) or one (5^) flux tubes is clearly 
shown in a*. 

Neither final state could be considered as a globally linear-force-free field but the individual 
twisted flux tubes do lie in regions of approximately constant a* . Thus while E^ is clearly in 
contradiction to a 'Taylor-like' relaxation (the Taylor state in that case being the uniform field), 
one could argue that the final state for S ^ has similarities with a Taylor state if we restrict ourselves 
to a domain over which significant current fragmentation occurs. For further considerations on 
the relevant Taylor state for 5 ^ we must consider the value of the total helicity for the field since 
this quantity is invoked as the constraint determining that state. 

Total helicity is used as a constraint under Taylor relaxation theory since, although it is not 
an invariant in a resistive MHD evolution, it is thought to be approximately conserved on the 
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relaxation timescales. To confirm this for these simulations we calculate the total relative helicity 
in the initial and final states for both fields. Here we use the reference field le^, which is the 
potential field satisfying the same boundary conditions. (Note that this is equivalent to calculating 
the total helicity in the torus obtained by identifying the top and bottom boundaries of our domain 
and assuming no flux links the hole of the torus.) In the initial states we find H(E^, t = 0) = 0.0 
and H{S^,t ^ 0) = 349.6. while for the final states we find H{E^, t - 350) = 1.2 and H{S^,t = 
350) = 350.9. Thus the conservation in both cases is excellent with only a small production of 
helicity (by magnetic reconnection) in each case. 

To see whether the final state for S ^ can be considered a Taylor state on the domain on which 
the current fragmentation occurs (see Figure |7]l we calculate the Taylor state on a cylindrical 
domain with the relevant parameters (radius r = 3, vertical flux of 9n and total helicity 350.8). 
The result is an axially symmetric Lundquist solution (Lundquist 1951) with constant a = 0.114 
so providing a reasonable qualitative and quantitative match to our results. 

Having described the basics of the resistive relaxations we proceed to consider the energetics 
of the process. This is of interest for determining whether and how relaxation processes such as 
these may heat the solar corona. 

4.2 Energy and heating 

In order to provide a solar coronal interpretation for these two relaxation events we present our re- 
sults in dimensional terms. To do so we choose characteristic values for the magnetic field strength 
as Bo = 10 G, the unit of length as /q - 1 Mm and the electron number density as 10^^ m^^. Ac- 
cordingly we have a unit of time as to = 1.45 s, temperature as To = 3.44 x 10^ K and velocity 
as vo - 689.7 km s"^ This gives our initial magnetic loops for both and a temperature of 
2.30 X 10^ K and a loop length of 48 Mm. Although the overall time of the relaxation then corre- 
sponds to 507 seconds, in considering this timescale it is important to note that the time taken for 
the full relaxation has previously (Pontin et al. 2011) been shown to be dependent on the value 
of the resistivity, with relaxation time increasing with decreasing resistivity. The resistivity taken 
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Figure 12: The total magnetic energy in excess of potential (upper image) and the total thermal 
energy above that of the initial state (lower image). In both cases the solid line is for while the 
dashed line is for 5^. 



as within reach of computing resources available to us remains several orders of magnitude too 
high when compared to assumed coronal values. Accordingly, the timescale for a real relaxation 
is likely to be greater than 500 s. 

The first thing that we note is that during the initial ideal relaxation process, the field 5 ^* is 
able to decrease its magnetic energy by more than . Specifically, both and begin with 
3.76 X lO'^/ of excess energy prior to the ideal relaxation but, after this ideal phase, E^ has 
1.25 X 10 7 remaining in excess of the potential field while 5^ has 1.06 x 10'^7. This is a result 
of the higher topological complexity of E^ - the highly braided field line mapping places a greater 
constraint on the ideal relaxation. Only once the resistive relaxation is begun, and the field begins 
to simplify, can this energy be liberated. 

The decay of magnetic energy in excess of potential during the resistive relaxation is shown 
in Figure 12 (upper image). While in the initial states for these resistive simulations E^ contained 
marginally more magnetic energy, over the course of the relaxation it releases 8.27 x 10^^7 of 
magnetic energy while releases only 4.85 x 10'^ 7. In terms of the fraction of free magnetic 
energy present at the start of each of the resistive simulations, 66.2% is released for E^ but only 
45.88% for . That is, E^ has been much more efficient at releasing energy. We claim this 
is because the higher degree of complexity in its initial state (as measured by the topological 
entropy) results in a much more fragmented, volume filling, system of current sheets (Figures [5} 
|6]and|7]l and so for a more efficient relaxation. 

Much of the magnetic energy released is converted into thermal form. The evolution of this 



quantity is shown in Figure 12 (lower image); we find an increase in thermal energy of 8.01 x 10^*^/ 
for E^ and 4.68 x 10'^7 for S^. The remaining energy takes the form of some residual non-zero 
kinetic energy in the final states of both systems. This is a result of large-scale weak oscillations 
in the flux tubes. 

We now examine changes in the temperature of both the E^ and S ^ loop systems during the 
relaxations. In doing so we must bear in mind that the initial plasma beta, p{t = 0) = 0.133, is 
relatively high and so temperature increases will be lower than if a lower plasma beta were used. 
Another factor is the lack of cooling terms in the energy equation which gives a competing effect 
in maintaining an artificially high temperature. Nevertheless we consider it useful to examine 
features of the temperature evolution in order to determine general trends. 

The initial temperature in both cases is 2.30 x lOf'K. In E^ the peak temperature increases to 
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Figure 13: Average temperature along field lines at their intersection with the z = plane (x,y,e 
[-3, 3]) for (upper panel) and (lower panel) at dimensionless times t = 40, 80, 120, 180 and 
350. 



2.76 X 10^ K and in 5^ to 2.92 x lO'^K. These are increases of factors 1.2 and 1.27 respectively. 
We have re -run the simulation for an initial plasma beta an order of magnitude lower {p{t = 
0) = 0.0133) and found an increase in maximum temperature of a factor 3.1 in that case. While 
the lack of realistic energy transport or loop stratification prevents us from making more realistic 
statements about heating along the loops, we can analyse the integrated temperature along field 
lines. This is justified by the excellent heat conduction along the magnetic field lines in the corona 
compared to the almost zero heat conduction perpendicular to the field. The quantity provides 
insight into the space-filling nature of the heating. The field line average temperature is shown 
in Figure 13 in the z = plane and for various stages of the evolution where the colour-scale for 
all images is normalised to the same maximum value (the end-state of S^). Some clear findings 
emerge here. The spatial distribution of high temperature is much more homogenous for E^. 
Nevertheless two brighter features are seen in the end-stages, corresponding to the two separated 
flux tubes in the end-state. In the early evolution several 'hot-spots' are present in the loop and 
these change location in time. We assert that the heating would be even more homogenous for 
lower resistivity as more current sheet fragmentation would be achieved (as documented in Pontin 
et al. 201 1). For heating primarily takes place in a thin layer in the centre of the domain with 
the temperature here being slightly higher than anywhere in the domain for E^ since the heating 
is less distributed. 

These contrasting behaviours in temperature distribution arise from a basic property of the 
initial configurations, namely, the difference in topological entropy. With higher topological en- 
tropy, E^ evolves through a highly fragmented system of current layers which allow for a uniform 
loop heating. The field S^, with lower topological entropy, evolves with fewer current layers 
which are more patchy in their distribution. The temperature profile is then correspondingly less 
homogeneous. 



5 Conclusions 

In this paper we have shown that the ability of non-ideal processes to heat coronal loops via mag- 
netic reconnection is crucially dependent on the nature of photospheric motions at the footpoint 
of the loops. Even for photospheric motions injecting the same amount of magnetic energy into 
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loop systems, the magnetic field configurations generated can lead to quite diff"erent amounts of 
energy being released in a relaxation event. 

Two eff'ects influence the amount and the distribution of heating. Firstly the heating is affected 
by the amount of energy that can be released in a relaxation on dynamic time-scales. Relaxation 
preserves the total helicity (Taylor, 1974) and possibly also further topological constraints (Yeates 
et al. 2010). Here, large amounts of magnetic helicity prove to be counter-productive as they raise 
the energy of the allowable end-state. In general more complex, non-coherent, braiding of the field 
lines allows for greater subsequent energy release. Secondly, the ability to release stored magnetic 
energy relies on a sufficiently fragmented, volume-filling system of current layers. Fragmentation 
is more efficient with more complex braiding, and also increases the homogeneity of the heating. 

f Note that the braid complexity (measured here by the topological entropy) is not entirely 
independent of the total helicity. High values of total helicity over a domain require a coherence 
in the structure that limits the topological entropy. Conversely, any non-trivial magnetic field on a 
closed or periodic domain must have a non- vanishing helicity density and so one can always find 
a sub-domain of such a field over which the total helicity is non-zero. Formalising a relationship 
between the two effects is beyond the scope of the present study. 

In order to support our claims we have presented simulations of the resistive MHD evolution 
of two magnetic loops. Both loops can be generated by pairs of rotational photospheric motions 
acting on footpoints of a uniform field and contain the same amount of magnetic energy. For our 
first field, E'^ , the two basic motions are in opposite directions while for our second field, , 
both motions are in the same direction. The boundary motions for lead to a field that is highly 
braided. The complexity in the field line mapping is visually apparent in a map of the squashing 
factor, Q (Figure [2]l and can be quantified by the topological entropy. The boundary motions 
for lead to a simpler pattern of field line connectivity with a simpler Q-structure (but similar 
maximum Q) and a lower topological entropy. 

The fields both undergo resistive relaxations on a fast, Alfvenic timescale. The relaxation 
for E^ is more efficient; more current sheets are generated which have a greater volume filling 
tendency. By the end of the relaxation the free magnetic energy of the field has been reduced 
by ~ 65% with this being converted primarily to thermal energy. The entire loop is heated in a 
relatively homogenous manner. The efficiency of the relaxation arises from the high topological 
entropy of the initial state. The relaxation for , with its lower topological entropy, is less 
efficient. A smaller number of current sheets result and these have a lower volume filling tendency. 
The free magnetic energy of the field is reduced by only ~ 45%. The heating of the loop takes 
place less uniformly, in one central location. 

Extrapolating to the solar corona, this mechanism provides a way to deposit heat in a spatially 
uniform way throughout the entire body of a coronal loop since complex braiding patterns of 
the field lead to a relaxation through a complex, space-filling set of current sheets. In order to 
determine how important this mechanism might be in practice, it would be useful to look at very 
high resolution data of surface motions (by way of fragment tracking) to consider the topological 
entropy of photospheric motions in different regions of the Sun over the solar cycle. Furthermore, 
a number of extensions to the models presented here should be made to allow for more detailed 
predictions relating to energetics. These include a reahstic stratification of the model atmosphere 
and additional physics to the energy equation. 
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